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Abstract 

We consider the motion of a two-dimensional body of arbitrary shape in a planar 
irrotational, incompressible fluid with a given amount of circulation around the body. 
We derive the equations of motion for this system by performing symplectic reduction 
with respect to the group of volume-preserving diffeomorphisms and obtain the rel- 
evant Poisson structures after a further Poisson reduction with respect to the group 
of translations and rotations. In this way, we recover the equations of motion given 
for this system by Chaplygin and Lamb, and we give a geometric interpretation for 
the Kutta-Zhukowski force as a curvature-related effect. In addition, we show that 
the motion of a rigid body with circulation can be understood as a geodesic flow on 
a central extension of the special Euclidian group SE(2), and we relate the cocycle in 
the description of this central extension to a certain curvature tensor. 

Contents 

1 Introduction 

2 Body-Fluid Interactions: Classical Formulation 

3 Body-Fluid Interactions: Geometric Approach 

3.1 Geometric Fluid-Solid Interactions 

3.2 Cotangent Bundle Reduction: Some General Results 



4 Derivation of the Chaplygin-Lamb Equations via Reduction by Stages 

4.1 Reduction by the Diffeomorphism Group 



1 Introduction 2 



4.2 Reduction by the Euclidian Symmetry Group 



18 



5 Geodesic Flow on the Oscillator Group 



20 



6 Conclusions and Outlook 



23 



A The Rigid Body Group 



24 



B The Group of Volume-Preserving Diffeomorphisms 



25 



1 Introduction 

We consider the motion of a rigid planar body, whose shape is not necessarily circular, im- 
mersed in a two-dimensional perfect fluid. We assume that the vorticity of the fluid vanishes, 
but we allow for a non-zero amount of circulation around the rigid body. This dynamical 
system is of fundamental importance in aerodynamics and was studied, among others, by 
Chaplygin, Kutta, Lamb, and Zhukowski; see Kozlov [1993], Kozlov [2003], and Borisov and 
Mamaev [2006] for a detailed overview of the literature and background information on this 
subject. A consequence of the non-zero circulation is the presence of a gyroscopic lift force 
acting on the rigid body, which is proportional to the circulation. This force is referred to 
as the Kutta- Zhukowski force. 

While various aspects of this dynamical system have been discussed throughout the fluid- 
dynamical literature (see for instance Milne-Thomson [1968] or Batchclor [1999] for a modern 
account), deeper insight into the Hamiltonian structure of this system had to wait until the 
work of Kozlov [1993] and Borisov and Mamaev [2006]. These authors identified a Hamil- 
tonian structure for the equations of motion found by Chaplygin and used this structure to 
shed further light onto the integrability of the system and to investigate chaoticity. This 
pioneering work raised several questions worth investigating in their own right. The first 
and foremost is that the non-canonical Hamiltonian structure of the rigid body with cir- 
culation is obtained by inspection, which immediately begs the question whether there are 
other, more fundamental reasons as to why this should be so. In this paper, we address this 
question following the approach of Arnold [1966]. In particular, we model the dynamics of 
the rigid body moving in a perfect fluid as a geodesic flow on an infinite-dimensional mani- 
fold. On this space, several symmetry groups act, and by performing successive reductions 
with respect to each of these groups, we obtain the Hamiltonian structure and equations of 
motion for the rigid body with circulation in a geometric way. This approach is also used 
in Vankerschaver, Kanso, and Marsdcn [2009] to derive the equations of motion for a rigid 
body interacting with point vortices. 

In this geometric approach, the symmetry groups of the fluid-solid problem are the particle 
relabeling symmetry group (the group of volume-preserving diffeomorphisms of the fluid 
reference space), and the special Euclidian group SE(2) of uniform solid-fluid translations 
and rotations. Symplectic reduction with respect to the former group is equivalent to 
the specification of the vorticity field of the fluid: in the problem considered here, this 
amounts to requiring that the external vorticity vanishes and that there is a given amount 
of circulation around the body. The group SE(2) acts on the system by combined solid- 
fluid transformations. After Poisson reduction with respect to this group, we obtain the 
Hamiltonian structure of Borisov and Mamaev [2006] . 

One straightforward advantage of the geometric description is that it provides immediate 
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insight into the dynamics, which might be harder to obtain by non-geometric means. As an 
illustration, we show that the symplectic leaves for the Poisson structure obtained by Borisov 
and Mamaev [2006] are orbits of a certain affine action of SE(2) onto se(2)*, which we 
compute explicitly. We show that these leaves are paraboloids of revolution, hence providing 
a geometric interpretation of a result of Chaplygin. 

An additional result of using geometric reduction is that we obtain a new interpretation 
for the classical Kutta-Zhukowski force. In particular, we show that the Kutta-Zhukowski 
force is proportional to the curvature of a natural fluid-dynamical connection, called the 
Neumann connection, which encodes the influence of the rigid body on the surrounding fluid. 
In this way, we exhibit interesting parallels between the dynamics of the rigid body with 
circulation and that of a charged particle in a magnetic field: both systems are acted upon 
by a gyroscopic force, the Kutta-Zhukowski force and the Lorentz force, respectively. The 
reduction procedure to recover the Hamiltonian description of Borisov and Mamaev [2006] 
then turns out to be similar to the way in which Sternberg [1977] derives the equations for 
a charged particle by including the magnetic field directly into the Poisson structure. 

Lastly, we also establish an analogue between the Kaluza-Klein construction for mag- 
netic particles and the dynamics of the rigid body with circulation. In the conventional 
Kaluza-Klein description, the dynamics of the magnetic particle is made into a geodesic 
motion by extending the configuration space and including the magnetic potential into the 
metric. Given the similarity between the Lorentz and the Kutta-Zhukowski force, a natu- 
ral question is then whether a similar geometric description exists for the rigid body with 
circulation. It turns out that this is indeed so: the extended configuration space in this 
case is a central extension of SE(2) by means of natural cocycle, which is again related to 
the curvature giving rise to the Kutta-Zhukowksi force. This extension of SE(2) is known 
as the oscillator group in quantum mechanics; see Streater [1967]. The rigid body with 
circulation is then described by geodesic motion on the oscillator group. 

Our description of the rigid body with circulation as a geodesic flow on a central extension is 
similar to the work of Ovsienko and Khesin [1987], who showed that the Korteweg-de Vries 
equation can be interpreted as a geodesic motion on the Virasoro group, an extension of the 
diffcomorphism group of the circle. It is also worth stressing that the similarity with the 
classical Kaluza-Klein picture cannot be taken too far: in the Kaluza-Klein description, one 
modifies the metric to take into account the magnetic field, whereas we leave the metric on 
SE(2) essentially unchanged and instead we deform the multiplication structure by means 
of a cocycle, giving rise to a central extension of SE(2). 

Outline of this Paper. We begin the paper by giving an overview of the classical litera- 
ture on fluid-structure interactions in section 2. In section 3 we recall some of the geometric 
concepts that arise in this context, most notably the particle relabeling symmetry group 
and the Neumann connection, and we give a brief summary of cotangent bundle reduction, 
which we use in section 4 to derive the Chaplygin-Lamb equations describing a rigid body 
with circulation. The geometry of the oscillator group and the link with the rigid body with 
circulation are explained in section 5, and we finish the paper in section 6 with a brief discus- 
sion of future work. In appendix A and B, some elementary results are collected about the 
geometry of the special Euclidian group SE(2) and the group Diff vo i of volume-preserving 
diff comorphisms . 
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2 Body-Fluid Interactions: Classical Formulation 

Consider a planar body moving in an infinitely large volume of an incompressible and inviscid 
fluid T at rest at infinity. The body B is assumed to occupy a simply connected region whose 
boundary can be conformally mapped to a unit circle, and it is considered to be uniform 
and neutrally-buoyant (the body weight is balanced by the force of buoyancy). Introduce 
an orthonormal inertial frame {ei^.3} where {ei,e2} span the plane of motion and e3 is 
the unit normal to this plane. The configuration of the submerged rigid body can then be 
described by a rotation 9 about e 3 and a translation x Q = x e\ + y e 2 of a point O (often 
chosen to coincide with the conformal center of the body) . The angular and translational 
velocities expressed relative to the inertial frame are of the form 9 and v — v x ei + v y e 2 
where v x — x , v y — y (the dot denotes derivative with respect to time t). It is convenient 
for the following development to introduce a moving frame {bi 2,3} attached to the body. 
The point transformation from the body to the inertial frame can be represented as 

„ „ „ fcosO — sin^A ,„ 

x = M + x , Re=[ sin9 cos0 ), (2.1) 

where x = x ei + y e 2 and X = Xhi + Yh 2 , while vectors transform as v = i?eV. The 
angular and translational velocities expressed in the body frame take the form fl — ilh^ 
(where f2 = 9) and V = V x h\ + V y h 2 (where V\ — x co&9 + y s'm9 and V 2 — —x sm9 + 
y cos9). Note that the orientation and position (9 7 x ,y ) form an element of SE(2), the 
group of rigid body motions in R 2 . The velocity in the body-frame £ = (ft, V x , V y ) T , where 
() T denotes the transpose operation, is an element of the vector space se(2) which is the 
space of infinitesimal rotations and translations in R 2 and is referred to as the Lie algebra 
of SE(2); for more details on the rigid body group and its Lie algebra, see Appendix A and 
references therein. 

Fluid Motion. Let the fluid fill the complement of the body in R 2 . The reference config- 
uration of the fluid will be denoted by J-o, and that of the body by Bq. The space taken by 
the fluid at a generic time t will be denoted by J- . Note however that as time progresses, the 
position of the body changes and hence so does its complement J 7 . In geometric mechanics, 
the mapping from the reference configuration J-o to the fluid domain T can be expressed 
as an element of the group of volume-preserving diffeomorphisms reviewed in Appendix B. 
In Section 3, we present an extension of this geometric approach to the solid- fluid prob- 
lem considered here but beforehand, we briefly describe, using the classical vector calculus 
approach, the fluid motion and the equations governing the motion of the submerged body. 

The fluid velocity u can be written using the Helmholtz-Hodge decomposition as follows 

u = V$ c + u v , (2.2) 

and has to satisfy the impermeability boundary condition on the boundary of the solid body. 
Now we explain the two terms in this decomposition. 

The potential function <E>£ is harmonic and represents the irrotational motion of the fluid 
generated by motion of the body. The subscript £ refers to the fact that $^ is determined 
by the body velocity £. 
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We emphasize that the motion of the body inside the fluid does not generate vorticity 
and only causes irrotational motions the fluid. The potential function $^ is a solution to 
Laplace's equation = 0, subject to the boundary conditions 



A$ c = 0, 



<9$ c 
dn 



(O x X + V) n, 



0. 



(2.3) 



<>B 



By linearity of Laplace's equation, one can write, following Kirchhoff (see Lamb [1945]), 

$ c = fi$ n + + V y %, (2.4) 

where <I>n, & x , <J? y are called velocity potentials and are solutions to Laplace's equation sub- 
ject to the boundary conditions on dB 



= (Xxn)-b 3 



OB 



dn 



= n • bi 



= n • b 2 



(2.5) 



OB 



OB 



The velocity u v is a divergence- free vector field and can be written as u v = V x \& + up, 
where V x \I> describes the fluid velocity due to ambient vorticity and Up describes the fluid 
velocity due to a net circulatory flow around the submerged body. The vector potential 
satisfies A\& = — u>, where uj = V x u v is the vorticity field, subject to the boundary 
conditions (V x $) ■ n = on dB and V x \P = at infinity. Clearly, in the absence of 
ambient vorticity, *f? is harmonic and can be written for planar flows as \& = ^e-s , where is 
referred to as the stream function and satisfies the boundary conditions (see Milne-Thomson 
[1968, §9.40]) 

®\ gB = V X Y - V y X - ^{X 2 + Y 2 ) = function of time only. (2.6) 



The harmonic vector field ur is non-zero only when there is a net circulatory flow around 
the body; it satisfies V-Up = and V x up = (i.e., Aup = 0) and the boundary conditions 
ur • n = on dB and up = at infinity. Note that, in three dimensional flows, one does 
not need the harmonic vector field ur- 1 

A harmonic stream function "fr associated with the circulation around the planar body can 
be defined such that ^rlgg = constant. The function ^ r can be found using a conformal 
transformation that relates the flow field in the region exterior to the body to that in the 
region exterior to the unit circle. For concreteness, let {X,Y) denote the body coordinates 
in the circle plane (which are measured relative to a frame attached to the center of the 
circle as a result of choosing the origin of the body frame in the physical plane to be placed 
at the conformal center). The stream function 'J'r can be readily obtained by observing 
that the effect of having a net circulation L around the body is equivalent to placing a point 
vortex of strength F at the center of mass of the body; namely, 

u r = V x (* r e 3 ), * r = -J- log (I 2 + Y 2 ). (2.7) 

47T 

Alternatively, the circulatory flow could be obtained as the gradient of a harmonic potential 
i»r (the harmonic conjugate to satisfying the Cauchy-Riemann relations). Note that $r 
would have to satisfy V$r ' n \os = 

1 In three dimensions, any closed curve in the exterior of a bounded body is contractible, so the harmonic 
vector field up may be set to zero. This result is due to the Poincare Lemma which can be alternatively 
stated as follows: a closed one-form on a (sub)-manifold with trivial first cohomology is globally exact. 
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Kinetic Energy. The kinetic energy of the fluid-solid system is simply the sum of the 
kinetic energies for both constitutive systems: 

1 f " l|2 jt/ , 1 ||Tr||2 , ! m 2 



T = T fluid + T body = - / Hull" + -m ||V|r + -m 2 . (2.8) 



where dV is a standard volume (more precisely, area) element on IR 2 . The kinetic energy 
of the rigid body can be readily rewritten in the form 

T body = ic T M fc C, M b := (I ^ . (2.9) 

where I is the 2-by-2 identity matrix. The kinetic energy of the fluid can be written as 

Tfluid 



If \\u\\ 2 dV=lf V($ c + $ r ) • V($ c + $ r ) dV 
Jjr Jjr (2.10) 

2 Jt jt 2 Jjr 

The first term in (2.10) can be rewritten using the divergence theorem, then employing (2.4) 
and (2.5), as follows 

\ jf V$ c ■ V$ c dV=\j^ c d ^dS= \C T M^ (2.11) 

where My is a 3 x 3 added mass matrix. Now, using the fact that and V$r are 

L2-orthogonal: since V$r = V x (f pes), where is uni- valued, we have 

f V$ c -V$ r rfV r =/ V$ c • V x (* r e 3 W = I $ c ^dS = 0, 
h h Job ds 

since \&r is constant on the boundary B, so that the tangential derivative d^r/ds vanishes. 
The last term in (2.10) can be treated as follows: in polar coordinates, &r{r, &) — r#, so 
that if we enclose the fluid-solid system in a large circular box of radius A, 

f A dr A 
V$ r • V$ r dV = 2irT 2 / — = 2tiT 2 log -. 
T Jr r R 

This constant term diverges logarithmically as A — > +oo. To remedy this, we regularize the 
kinetic energy by discarding this infinite contribution - a remedy also used in Lamb [1945]; 
Borisov and Mamaev [2006]. That is, we consider the kinetic energy of the solid-fluid system 
to be given by 

T=ic T (M 6 +M / )C. (2.12) 



Equations of Motions. The equations governing the motion of the body in potential 
flow with non-zero circulation around the body but in the absence of ambient vorticity are 
of the form 

rr = (P x v) • b 3 

. V ' (2.13) 
P = PxO + rb 3 xV 

Here, II and P denote the angular and linear momenta of the solid-fluid system expressed 
in the body frame. They are given in terms of the velocity in body frame by (here, T is 
given by (2.12)) 

dT dT 

" = " .v < 214 > 
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One of the main objectives of this paper is to use the methods of geometric mechanics 
(particularly the reduction by stages approach) to derive equations (2.13) governing the 
motion of the body in potential flow and with non-zero circulation. The case of a body of 
arbitrary geometry interacting with external point vortices is addressed in Vankerschaver, 
Kanso, and Marsdcn [2009]. It would be of interest to extend these results to the case of a 
rigid body with circulation moving in the field of point vortices. 

3 Body-Fluid Interactions: Geometric Approach 

In this section, we first establish the structure of the fluid-solid configuration space as a 
principal fiber bundle and show that there exists a distinguished connection on this bundle. 
We then recall some general results from cotangent bundle reduction that will be useful to 
reduce the system and derive the equations of motion (2.13). 

3.1 Geometric Fluid-Solid Interactions 

The Configuration Space. We describe the configurations of the body-fluid system by 
means of pairs (<?,</?), where g is an element of SE(2) describing the body motion and 
ip : — > R 2 is an embedding of the fluid reference space Fq into R 2 describing the fluid. 
The pairs (g, ip) have to satisfy the impermeability boundary conditions dictated by the 
inviscid fluid model. The embedding ip represents the configuration of an incompressible 
fluid and has therefore to be volume-preserving, i.e., ip*(dV) = dVo, where dVo and dV are 
volume forms on T and K 2 , respectively. We denote the space of all such volume-preserving 
embeddings by Emb vo i(J"o, ^ 2 ). We denote by Q the space of all pairs (g,(p), g £ SE(2) 
and <p £ Emb vo i(J*o,K 2 ) that satisfy the appropriate boundary conditions. That is, the 
configuration manifold Q of the body-fluid system is a submanifold of the product space 
SE(2) x Emb vol (Jo,IR 2 ). 

Tangent and Cotangent Spaces. At each (g,<p) £ Q, the tangent space Tr g ^)Q is a 
subspace of T g SE(2) x T 1/3 Emb vo i(.7 : o,M 2 ) whose elements we denote by (g, <p,g, <p). Here, g 
is an element of T g SE(2) and <p is a map from To to TR 2 such that (p(x) £ T V ^R 2 for all 
x £ Fq. Note that g represents the angular and linear velocity of the rigid body relative to 
the inertial frame while <p represents the material or Lagrangian velocity of the fluid. It 
is easier however to represent the elements of TQ using the rigid body velocity expressed in 
body frame ( and the fluid velocity u. Note that, in the group theoretic notation, the body 
velocity may be defined as ( = g~ x g and the spatial or Eulerian velocity field of the 
fluid may be defined as u :— ip o p~ x . The vector field u is a vector field on T := p(J-"o), in 
contrast to (p, which is merely a map from to TR 2 . We emphasize that £ and u cannot 
be chosen arbitrarily, but have to satisfy the impermeability boundary conditions. 

The cotangent space T* g at a point (g,<p) G Q consists of elements (g,(p,n,a), where 
7r = MhC £ se(2)* is the momentum of the submerged body and a £ ^ 1 (^ r ) is the one-form 
dual of the velocity field u, see Appendices A and B for more details. 

Kinetic Energy on T*Q . The kinetic energy in (2.12) can be used to define a metric 
on the cotangent bundle T*Q. To verify this, it is informative to recall here the Hodge 
decomposition of differential forms (see Abraham, Marsdcn, and Ratiu [1988]). Any one- 
form a £ ^{J 7 ) can be decomposed in a unique way as a = + ap, where $ £ r2°(J r ), 
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W € il 2 (J r ), and «r is a harmonic form: dar = <5c*r = 0. The Hamiltonian of the solid-fluid 
system can now be defined as the function H : T*Q M given by 

H(g, <p, 7r b , a) = \ ||a|| 2 + \ix T WT^ tt = \ ||d$ c || 2 + i ||<Sv]/|| 2 + ^ Mr 1 7r fe , 

where the norm in ||a|| 2 is that in (B.f) on one-forms induced by the Euclidian metric 
and where we have again discarded the infinite constant ||ar|| • Here we have denoted the 
momentum by 7Tb to emphasize the fact that this variable encodes the momentum of the rigid 
body only. Later on, we will consider the momentum of the combined solid-fluid system, 
which will be denoted by 7r = (II, P). 

The first term in the right-hand side of the Hamiltonian can be expressed in terms of the 
added mass matrix: 

||d$ ? || 2 = / d<£>£ A *d$£ = (t> * d<&c - / * d^c 
Jr JdB Jr 

= 1 <P c ^dS- [ $ c V 2 $ c dV = C T M f C, 

JdB dn Jr 

so that the Hamiltonian becomes 

H(g,ip,7r b ,a) - l - ||(^|| 2 + ^M^Mfc + M/)Mr x tt 6 (3.1) 

using Q = M^ 1 ^. The term involving 5^ yields the kinetic energy due to vortical structures 
present in the fluid and is zero for the rigid body with circulation. 



The Action of the Group of Volume-Preserving Diffeomorphisms. The group of 
all volume-preserving diffeomorphisms T)\Q vo \{F) acts from the right on Q by composition: 
for any (g, tp) € Q and 4> G Diffvol^) we define 

(g,<p) ■ <t> ■= (g,tp°<P)- 

This action leaves the kinetic energy (3.1) on Q invariant, since the Eulerian velocity field u 
is itself invariant. The Diff vo i(J r )-invariance represents the particle relabeling symmetry . 
The manifold Q is hence the total space of a principal fiber bundle with structure group 
Diff vol (.F) over SE(2). Here, the bundle projection projection pr : Q — ¥ SE{2) is simply the 
projection onto the first factor: pr(g, tp) — g. One can readily show that the infinitesimal 
generator Xq corresponding to an element X S Xvoi^o) is given by 

X Q (g,^):=(0,^X)€T M Q. (3.2) 



The Momentum Map of the Particle Relabeling Symmetry. The group Diff vo i(.Fo) 
acts on Q and hence on T* Q by the cotangent lifted action. We now compute the momentum 
map corresponding to this action, i.e. the particle relabeling symmetry. This is a map J 
from T*Q to 3£* ol (J" ), and we recall from appendix B that X* ol (J" ) = dfi 1 (J r ) x M. 
Consequently, the momentum map has two components, corresponding with circulation and 
vorticity (pulled back to the reference configuration To). The statement that the momentum 
map is conserved then translates into Kelvin's theorem that the vorticity is advected with 
the fluid, and that the circulation around each material loop is conserved. 

Proposition 3.1. The momentum map J of the DiS vo \(J-o)- action on T*Q is given by 



J(TVb,a) = (d</?*a,r), where T = 
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Proof. We use the well-known formula for the momentum map of a cotangent lifted action 
(see Abraham and Marsden [1978]). For each X g X vo \(To), we have 

(J(ir b ,a),X) = ((ir b ,a),X Q (g, l p))= [ (a, <p.X) dV = f (<p* a, X) dV , 

so that J(nf,,a) — [tp*a] € £* ol (.Fo). After composition with the isomorphism (B.4) we 
obtain the desired form (3.3). ■ 

The One- Form ar with Circulation T. Recall from section 2 that having circulation T 
around the rigid body is equivalent to placing a point vortex of strength T at the conformal 
center of the body, whose velocity field up was given in (2.7). 

In the remainder of this paper it will be easier to work with the one-form ar on T given by 
ar = Up, or explicitly by 

a r = 8(V T dV), (3.4) 
where is the stream function (2.7), and dV is the volume form on R 2 . 
Proposition 3.2. The one-form ar is a harmonic one-form on J- and satisfies 

I a r = r. (3.5) 

JOB 

In particular, ar is Li-orthogonal to the space of exact one-forms. 

Proof. The form ar satisfies Sar = by definition, and we have that dar = TS(X.)dV, so 
that dar = in the fluid domain J- '. Hence, ar is harmonic and (by means of the Hodge 
theorem) L2-orthogonal to the space of exact one-forms. The line integral (3.5) follows from 
the expression (2.7) for the stream function $r- B 

The Action of the Euclidian Symmetry Group. In addition to the right principal 
action of Diff vo i(.7-o) on Q described before, the special Euclidian group SE(2) acts on Q by 
bundle automorphisms from the left. In order words, there is an action ip : SE(2) x Q — >• Q 
given by 

(g, <p))=h - (g, (f) = (hg, hip), (3.6) 

for all h € SE(2) and (g, <p) € Q. The embedding h<p is defined by (hip)(x) = h ■ <p(x), 
where the action on the right-hand side is just the standard action of SE(2) on M. 2 . From 
a physical point of view, the S'i?(2)-action corresponds to the invariance of the combined 
solid-fluid system under arbitrary rotations and translations. 

The Neumann Connection. The bundle pr : Q — >• SE(2) (defined by the Diffvo^J 7 )- 
invariance) is equipped with a principal fiber bundle connection, which was termed the 
Neumann connection in Vankcrschaver, Kanso, and Marsden [2009]. This connection 
seemed to have appeared first in Lewis, Marsden, Montgomery, and Ratiu [1986] (see also 
Koiller [1987]) and essentially encodes the effects of the rigid body on the ambient fluid. 

The connection one- form A : TQ — > 3t vo i (Fo) of the Neumann connection is defined in terms 
of the Helmholtz-Hodge decomposition (2.2) of vector fields: if (g, <p, £, u) is an element of 
T(g lV> )Q, then 

•A(g, <p, C, u) = (/3*u V) (3.7) 
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where u v is the divergence-free part of the Eulerian velocity u in the Helmholtz-Hodge 
decomposition (2.2). It can be shown (see Vankcrschavcr, Kanso, and Marsden [2009]) that 
A satisfies the requirements of a connection one-form, and that A is invariant under the 
action of SE(2) on Q: 

Aj#)( T ^(fi u )) = As,v)(C' u ) 

for all (C,u) e T( 3iV )Q. Here iph = ^(h, ■), with tp the S*i?(2)-action (3.6). Given the exact 
form of the Neumann connection, one can compute its curvature, which is a two-form B 
on the total space Q with values in JEvoip^o)- It turns out that there exists a closed- form 
formula for the curvature, which was first determined by Montgomery [1986] and further 
generalized by Vankerschaver, Kanso, and Marsden [2009]. More precisely, we compute an 
expression for the contraction (//, B), where \x is an arbitrary element of the dual space 

Proposition 3.3. Let (d,Ui) and (C 2 ,u 2 ) be elements ofT^ gtp )Q and denote the solutions 
of (2.3) associated to Ci resp. £ 2 by $i an d 3> 2 . Let /i be an arbitrary element of X* o1 (Fq). 
Then the fi-component of the curvature B is given by 

(M,B( g ^)((Ci,u 1 ),(C 2 ,u 2 ))> = « M ,d$ 1 Ad$ 2 »- I aA*(d$!Ad$ 2 ), (3.8) 

JOB 

where ((•, •)) is the metric on the space of forms on T defined in (B.l). 

In what follows, it will be necessary to have an expression for the curvature in terms of 
the elementary stream functions, rather than the elementary velocity potentials. Such an 
expression can be easily obtained by noting that the stream function f is a harmonic 
conjugate to the velocity potential In particular, if $i,<fr 2 are the velocity potentials 
introduced in the statement of proposition 3.3, with their associated stream functions ^t 2 , 
then d^! A d , 5 2 = d$ x A d$ 2 . The /Lt-component of the curvature (3.8) hence becomes 

{H,B) = «/i,d*i Ad* 2 )) - f a A *(d*i Ad* 2 )- (3.9) 

JOB 

We established the structure of the fluid-solid configuration space Q as the total space of 
a principal fiber bundle with structure group Diff vo i(.7 : o). We also showed that there exists 
a distinguished connection on this bundle, which is invariant under the action of SE(2) on 
Q. In order to reduce the system and derive the equations of motion (2.13), we will need 
the framework of cotangent bundle reduction, which we now describe. More information 
and proofs of the results quoted below can be found in Marsden and Perlmutter [2000] and 
Marsden, Misiolek, Ortega, Perlmutter, and Ratiu [2007]. 

3.2 Cotangent Bundle Reduction: Some General Results 

We describe cotangent bundle reduction in a general context. Let the unreduced configura- 
tion space be denoted by Q and assume that a Lie group G acts freely and properly on Q 
from the right, so that the quotient space Q/G is a manifold. Furthermore, we assume that 
the quotient space Q/G is equal to a second Lie group H, which acts on Q from the left. In 
Section 4, Q will be the fluid-solid configuration space, G will be the group Diffvoi^o) and 
H will be the special Euclidian group SE(2). Our first goal is to reduce by the G-action 
and describe the reduced phase space. This is described in theorem 3.4. In the second stage 
of the reduction, we will then do Poisson reduction with respect to the residual group H in 
order to obtain a Poisson structure on the twice reduced space. This is described below in 
theorem 3.5. 
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Recall that the quotient projection prg G : Q — > Q/G defines a right principal fiber bundle, 
and assume that a connection on this fiber bundle is given, with connection one-form A : 
TQ — > q. For more information about principal fiber bundles and connections, see Kobayashi 
and Nomizu [1963]. The group G acts on T*Q by cotangent lifts, and we denote the 
momentum map of this action by J : T*Q — > g*. The next theorem characterizes the 
reduced phase space (T*Q) fl := J _1 (/i)/G M , where fi is an arbitrary element of g*. Here, 
G p is the isotropy subgroup of /i, defined as follows: g € G M if Ad*fi = /i, where Ad* 
denotes the coadjoint action of Gonj*. The reduced phase space (T*<5) M can be described 
in full generality, but we focus here on the special case where the isotropy group G M is the 
full symmetry group: G M = G. We will see in section 4.1 that this is the relevant case to 
consider for the rigid body with circulation. 

Theorem 3.4. Let G be a group acting freely and properly from the right on a manifold 
Q so that prg G : Q — > Q/G is a principal fiber bundle. Let A : TQ — > g be a connection 
one-form on this bundle. Let ft£g* and assume that G M = G. 

Then there is a symplectic diffeomorphism between {T*Q)^ and T*(Q/G), the latter with 
symplectic form w can — S„; here w can is the canonical symplectic form on T*(Q/G) and 
= prg^ G /3 M , where Wq/g '■ T* {Q/G) — >• Q/G is the cotangent bundle projection, and /3 M 
is determined through 

w ^ G ^ = d(ti,A). (3.10) 

Outline of the Proof. This is a special case of theorem 2.3.3 in Marsden, Misiolek, Or- 
tega, Pcrlmuttcr, and Ratiu [2007]. We just recall the explicit form of the isomorphism 
between (T'*Q) M and T*(Q/G); the proof that this map also preserves the relevant symplec- 
tic structures can be found in Marsden, Misiolek, Ortega, Perlmuttcr, and Ratiu [2007]. 

The isomorphism ip^ : (T*Q) fl -> T*(Q/G) is the composition of the map shifty : {T*Q)^ — s- 
(T*Q) and the map ip : (T*Q) Q ->■ T*(Q/G): 

= Vo ° shifty. (3.11) 

Both of these constitutive maps are isomorphisms. The map shifty is defined as follows: we 
first introduce a map Shifty : J~ 1 (fj,) — > J _1 (0) by 

Shift AI (a 9 ) = a q - (u,A(q)) ■ 

It can easily be verified that Shift is G-invariant so that it drops to a quotient map shifty : 

(r*Q) M ->• {T*Q)o- 

Secondly, the map ip : (T*Q) -t T*(Q/G) is defined by noting that 

J-^O) = {a q £ T*Q : (a q ,Z Q (q)) - for all £ € fl} 

so that the map tp : J _1 (0) -> T*(Q/G) given by 

(ip {a q ),Tn Q ^ G (v q )) = (a q ,v q ) (3.12) 

is well defined. The map <po is easily seen to be G-invariant and surjective, and hence 
induces a quotient map ip : {T*Q) T*(Q/G). ■ 

Note that the isomorphism between (T*Q) M and T*Q/G is connection-dependent. As a 
result, the reduced symplectic form on T*Q/G is modified by the two- form /3 M , which is 
traditionally referred to as a magnetic term since it also appears in the description of a 
charged particle in a magnetic field (see Guillemin and Sternberg [1984]). 
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Having described the reduced phase space (T*Q) M , we now take into account the assumption 
made earlier that the base space Q/G has the structure of a second Lie group H , which acts 
on Q from the left and leaves the connection one-form A invariant. In this case, the reduced 
phase space (T*Q)^ is equal to T*H, equipped with the magnetic symplectic structure 
described before. As a result, H acts on T* H, and it can be checked that this action leaves 
the symplectic structure invariant. It would now be possible to do symplectic reduction as 
above for the if-action as well to obtain a fully reduced symplectic structure. However, all 
that is needed in section 4.2 is an expression for the reduced Poisson structure, described 
in the following theorem. 

Theorem 3.5 (Theorem 7.2.1 in Marsden, Misiolck, Ortega, Perlmutter, and Ratiu [2007]). 
The Poisson reduced space for the left cotangent lifted action of H on (T*H, w can — B^) is 
t)* with Poisson bracket given by 

{/,3}b(m) = - (ft, 

for f,geC°°(t)*)- 

The theorem in Marsden, Misiolek, Ortega, Perlmutter, and Ratiu [2007] is proved for right 
actions, whereas the action of H here is assumed to be from the left. However, the same 
proof continues to hold, mutatis mutandis. 

Lastly, we recall the definition of the £?f)-potential (see Marsden, Misiolek, Ortega, Perlmut- 
ter, and Ratiu [2007]), which plays the role of momentum map relative to the magnetic term. 
While this potential can be defined on an arbitrary manifold with a magnetic symplectic 
form, we treat just the case of a Lie group H with left-invariant symplectic form u> can — B^ 
as before. 

Definition 3.6. Let H be a Lie group with left-invariant symplectic form w can — B^ . Suppose 
there exists a smooth map tp : H — » t)* such that 

i 6B B ll = d(^,0, (3-14) 

for all £ £ f). Then the map ip is called the B\)-potential of the H -action, relative to the 
magnetic term B^ . 

In what follows, we will always assume that tp exists. In this case, ip is defined up to an 
arbitrary constant, which we normalize by assuming that tp(e) = 0. Under this assumption, 
the non-equivariance one-cocycle a : H — > f) associated to tp, 

*(g) = ip(g)-Ad* g -^(e), (3.15) 

coincides with tp. We will make no further distinction between ip and a. The importance 
of ip lies in the fact that we may characterize the symplectic leaves of the magnetic bracket 
(3.13) in t)* as orbits of a certain affine action of H on f)*. 

Proposition 3.7. Let (T*H,uj can — B^) be a He group with a magnetic symplectic form. 
The symplectic leaves in f)* of the magnetic Poisson structure (3.13) are the orbits of the 
following affine action: 

g-H = Ad* g -ifi + ip(g), (3.16) 

where g £ H, [i 6 \f . 
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Proof. See theorem 7.2.2 in Marsden, Misiolek, Ortega, Perlmutter, and Ratiu [2007]. ■ 



4 Derivation of the Chaplygin-Lamb Equations via Reduction by Stages 



13 



4 Derivation of the Chaplygin-Lamb Equations via Re- 
duction by Stages 

4.1 Reduction by the Diffeomorphism Group 

We impose the condition that the fluid has constant circulation T by considering the sym- 
plectic reduced space J^ 1 (0, r)/Diff vo i(J r ). Our goal is to use cotangent bundle reduction 
to express this space in a more manageable form (via establishing an isomorphism that is 
connection-dependent) and to obtain an explicit expression for the reduced symplectic form 
on this space. 

The reduction theorem 3.4 deals with the case where the isotropy group coincides with the 
full group. It can be easily verified that this condition is satisfied for the rigid body with 
circulation: here, the relevant momentum value is (0,r), and by (B.5) we have that for all 
e DiflVoiGFo), 

CoAd (o,r) = (o,r), 

so that (Diff vol (J r o))(o,r) = Diff vol (J r )- 

The Reduced Phase Space. As an introduction to the methods of this section, we 
establish an isomorphism <pr between the reduced phase space J _1 (0, r)/Diff vo i(J r o) and 
T*SE(2) in a relatively ad-hoc manner. We then show that cotangent bundle reduction 
yields precisely this isomorphism. 

Let us first introduce the linear isomorphism m : se(2)* — > se(2)* given by 

m(7r b ) = nb + MfM^TTb = (M 6 + M / )M^ 1 7r 6 (4.1) 

where and My are the body mass matrix and the added mass matrix, respectively. The 
map m will prove to be crucial later on; its effect is to redefine the momentum of the rigid 
body in order to take into account the added mass effects. 

We now derive an explicit expression for the level sets J _1 (0) and J _1 (0, T) of the vorticity 
momentum map. Let (g,(p,Wb,a) be an element of T*Q: the requirement that J(iTb,a) — 
(0, r) is equivalent to 

da = and <p ip* a = Y. 
JdB 

By means of the Hodge decomposition, this implies that 

a = or + d$£, 

where ar is given by (3.4) and $^ is the solution of the Neumann problem (2.3) with 
boundary data £ = M^ 1 ^. Hence, we may define an isomorphism tpr from J~ 1 (0, T) to 
Q x se(2)*, given by 

tpr ■ {g,<p,irb,a) ^ (g,f,^(^b)) € Q x st(2)*, 

where m is the map (4.1). Likewise, there is an isomorphism i/jq : J _1 (0) — > Q x se(2)*, 
obtained by realizing that J _1 (0) consists of elements of the form (g, tp, %bi where $^ 

has the same interpretation as above. 

The map ipr is easily seen to be Diff vo i(J r o)-equivariant and hence drops to a quotient 
isomorphism <pr between J _1 (0, r)/Diff vo i(J : o) and T*SE(2), given by 

<pt ■ [(g, n, ")] >->• (jg,tn(7r b )). 

We see that the effect of the map ip^ is to eliminate the influence of the fluid entirely, except 
for the added mass effects, which are encoded by m. 
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The r-Component of the Connection One-Form. In order to compute the magnetic 
term (3.10), we need an expression for the connection one-form, contracted with the momen- 
tum value at which we do reduction. For the rigid body with circulation, the connection 
is the Neumann connection, and the momentum value is (0,T). Therefore, let A be the 
connection one-form of the Neumann connection as in (3.7) and define the T-component of 
A to be the one-form At ■ TQ —> R given by 

A r (g, V, C, u) := ((0, T),A(g, if, C, u)) , 

where on the right-hand side we interpret (0, T) as an element of % vo i(J-q)* . In other words, 
Ar(g, <p, C) u ) computes the divergence-free part of u and contracts it with the element (0, T) 
of X vo i(J-"o)*- This hints at the fact that At is nothing but the one-form ar of (3.4), which 
we prove in the next proposition. 

Proposition 4.1. The T-component of the connection one-form A is related to the form 
ar by the following relation: for all (g,(p) € Q and (£, u) £ T( g ^)Q, we have that 

■4 r (C,u) = J^a T {u)dV. (4.2) 

Proof. By definition, we then have that At is given by 

^(C,u) = ((0,r),^(C,u)> = (a r ,u v ) = J or(u) dV, 

where we have used the Helmholtz- Hodge decomposition u = u v + , together with the 
fact that ar satisfies (3.5) and is ^-orthogonal to gradient vector fields. ■ 

The T-Component of the Curvature. As a second step towards the computation of the 
magnetic symplectic form (3.10), we need a convenient expression for the exterior derivative 
d„4r- While it is in theory possible to compute the derivative directly from the expression 
(4.2) for Ar, we can avoid much of the technicalities by means of the following insight: for 
the rigid body with circulation, d.4r is nothing but the T-component of the curvature of 
the Neumann connection, as defined below. Recall that an explicit expression for the latter 
was established in proposition 3.3. We define the T-component of the curvature B as before 
by the following prescription: 

Br((C,u),(^v)):=((0,r) J B((C,u),(C,v))) J 

again relying on the fact that (0, T) e X V oi(Fo)* , while B(((, u), (£, v)) is an element of 

Proposition 4.2. The T-component of the curvature is a left SE (2) -invariant two-form on 
SE(2) given at the identity by 

B r (e)=Te* x Ae* y . (4.3) 

Proof. The left-invariance of Br follows from the fact that the Neumann connection itself 
is left invariant under the action of SE(2) upon itself. We may hence restrict our attention 
to the value of the curvature at a point (e, tp) G Q, where e is the identity in SE(2). The 
fact that Br drops to SE(2) will be obvious once we determine its exact form, but can also 
be proved directly; see for instance Vankcrschaver, Kanso, and Marsden [2009]. 

The first integral in the expression (3.9) for the curvature vanishes since the integration is 
over the fluid domain, where dar = 0. The second integral can be computed explicitly, 
since we only need the expression (2.6) for the stream function on the boundary. 
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Let (g,ip,Q,Ui), i — 1,2, be elements of T( g , v )Q and write the velocity as d = (f2;,Vj), 
where V, = C/j(cosaj,sinaj). Consider the stream functions given in (2.6) corresponding 
to the rigid body motions (fti, V^). On the boundary dB, we then have that 

* (d^i A d<3>2) = U1U2 sin(o;i — 02) + fiiL^zcosa^ + y sin 02) 

— ^2U\{x cosai + ysinai), 

so that the curvature is given by 

B r ((Ci, Ui), (C2, u 2 )) = - / a>r A *(d*i A d* 2 ) = -TU^ sm(a!i - a 2 ). 

JOB 

Here, we have used the fact that 

ar — r and d> xar = f yotr = 0. 
as Job Job 

Note that the value of Br does not depend on ui, u 2 so that Br drops to SE{2) as claimed. 
Finally, using the basis (A. 6) of se(2)*, we may write the curvature at the identity as 

6 r (e) = re:Ae;. 

This concludes the proof. ■ 

We now prove the main result of this section, that the exterior derivative d„4r is equal to the 
T-component Br of the curvature. More precisely, recalling the projection pr : Q — > SE(2), 
we have 

pr*B r = cL4 r . (4.4) 

Note that this is not true for arbitrary reduced cotangent bundles, and that this is highly 
specific to the case of rigid bodies with circulation. The expression (4.4) can be proved as 
follows: because of the Cartan structure formula B = d.4+ [.4, *4] for right actions, we have 
that 

d^l r =pr*B r -<(0,r) ! [A^]>. (4.5) 

It remains to show that the second term on the right-hand side vanishes. It can be shown 
(see the proof of theorem 4.2 in Marsden and Wcinstcin [1983]) that, for divergence-free 
vector fields u Vi i and u v , 2 tangent to dF and arbitrary one-forms a, the following holds: 

J a([u v ,i, u v , 2 ]) dV = J da(u v ,i, u V)2 ) dV. (4.6) 

If we now put u Vj i := A(Q, Uj), i = 1,2, we may rewrite the second term in (4.5) as follows: 
((0,r),L4(Ci,ui),.A(C2,u 2 )]) = ((0,r) l [u v ,i > u v ,a]) 

= J a r ([uv,i,u v ,2])^ = J dar(u v ,i,u V)2 )dV, 

where ar is the one-form given by (3.4). However, the integral on the right-hand side is 
again zero because the integration is over J- ' , while the support of dar is concentrated 
inside the rigid body. This shows that d„4r in (4.5) is nothing but the T-component of the 
curvature of the Neumann connection, as calculated in proposition 4.2. 
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The Reduced Phase Space. We formally establish the isomorphism tp? of theorem 3.4 
between the reduced configuration space J _1 (0, r)/Diff vo i(J r ) and T*SE(2). Recall that 
(pr is defined as the composition cpo o shiftr , where tpo and shiftr were defined in the proof 
of theorem 3.4. The following two lemmas are devoted to the explicit form of these two 
constitutive maps for the rigid body with circulation. 

Lemma 4.3. Consider the map m : se(2)* — > se(2)* defined in (4-1)- The map (po : 
(T*Q)q — > T* SE(2) defined in the proof of theorem 3.4 is given in terms of m by the 
following expression: 

tp {{g,ip,n b ,d<S> c ]) = (g,m(TT b )). (4.7) 

Proof. Consider an element (g, ip, d$^) of J _1 (0). In a left trivialization, the map 
<p : J _1 (0) T*SE{2) defined in (3.12) is given by 

(£(7T 6 , d* c ), = ((7Tb, d$ c ), (£, V$ c )) = (tt 6) + (d$ c , V0> c ) , 

where now $j is the solution of (2.3) with boundary data £. The right-hand side can then 
be written as follows, using the definition of the added-mass matrices: 

((7r b ,d$ c ),(£,V$ e )) = ( 7 r 6 ,0+^d$ c .V$ ? dF = C T (M b +M / )e 
= ((I + M/M^ 1 )^) 
since 7r b = M b £. In other words, we have that 

(<p(g, ip, 7T b , d$ c ), = (m(7T b ),0 , 
so that ip, n b , d$ c ) = (3, m(7r&)). ■ 

We now determine the map shiftr : (^*Q)(o,r) ~> {T*Q)o- Because of proposition 4.1, the 
map Shiftr : J~ 1 (0,T) — > J -1 (0) is in our case a simple shift by ctr- 

Shift r (g, <p,7r b ,a) = {g,(p,ir b ,a- a r ). 

Physically speaking, we may think of a £ J _1 (0,r) as a one- form with circulation T; by 
subtracting ar from a, we obtain a one-form with zero circulation. For the quotient map 
shiftr : (T*Q)( 0j r) — > (T*Q) we then have the following result. 

Lemma 4.4. The quotient map shiftr : (T*Q)m t r) —> (T*Q)q is given by 
shiftr([(ff, 95,7Tb, a)]) = [(3, 93, ?r b , a - ar)]. 

By concatenating these two maps, it is now straightforward to find the explicit form of the 
isomorphism tpr between (T*Q)(o,r) an d T*SE(2): from (3.11) we have that the isomor- 
phism is given by 

ipr : [(g, <p, 7r b , a)} i-» {g, m(n b )). (4.8) 

The general theory of cotangent bundle reduction guarantees that this map is an iso- 
morphism, but it is instructive to construct the inverse mapping explicitly. For every 
(<7, 7r) <E T*SE(2), put 7r b := m _1 (7r) and set a := ar + d$^, where is the solution 
of the Neumann problem (2.3) with boundary data £ = M^" 1 ^. Furthermore, choose an 
arbitary fluid embedding ip such that (g, ip) 6 Q. The inverse mapping ip^. 1 is then given by 

ip^ 1 (g,ir) = [(g,ip,n b ,a)], where n = m _1 (7r b ) and a = ar + d<I>^. (4-9) 

It is straightforward to check that ip^ 1 o tpr = ipr o ip^. 1 = id. 
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The Reduced Hamiltonian. As a last step towards establishing a reduced Hamiltonian 
formulation for the rigid body with circulation, we need to find an appropriate expression 
for the Hamiltonian on this space. This can be done by computing first the Hamiltonian on 
the unreduced space from the kinetic energy (2.8) and then noting that the result induces 
a well-defined function on the reduced phase space. 

Proposition 4.5. The Hamiltonian on the reduced phase space T*SE(2) is given by 

H red ( 5 ,7r) = i 7 r r M- 1 7r, (4.10) 

where M = M b +M f . 

Proof. The proof relies on the explicit form of the isomorphism ipr in (4.8). Let (g, ir) be 
an element of T* SE(2) and consider an arbitrary element (g, ip, Kb, ot) £ J~ (0, T) such that 
Vr 7r ) = 1(5' fi nb ' a )J- R- ecan from the discussion following (4.9) that 7T(, and a are given 
by 

TTb = m 1 (tt) and a = ar + d$£, (4-11) 

where is the solution of the Neumann problem (2.3) with boundary data £ = M^" 1 ^. 
The relation between the reduced Hamiltonian -ff rc d on T* SE(2) and the Hamiltonian (3.1) 
is then written as 

H Ted (g,Tr)=H(g,<p,Tr b ,a) = -tr^ M b _1 (M; ) + M^M^ 1 7T(,, 

where we have used the fact that the fluid is irrotational, so that U 7 = in (3.1). Keeping 
in mind that tt = m(7Tf,), or alternatively that n = (Mb + My^M^ -1 ^, we finally obtain the 
following expression for the reduced Hamiltonian: 

Hrcd{g,n) - l -Tt T (Mfc+M/)- 1 ^, 

which is precisely (4.10). ■ 



Summary. In this section, we used the framework of cotangent bundle reduction to obtain 
expressions for the reduced phase space, its symplectic structure, and Hamiltonian after 
reducing with respect to the particle relabeling symmetry. We summarize these results in 
the theorem below. 

Theorem 4.6. Let rel and consider the associated element (0,T) 6 3E* ol (.Fo). Then the 
following properties hold: 

• The isotropy group Diff vo i(J r o)(o i r) o/(0,T) is the whole of the group Diff vo i(.7 : o). 

• The reduced phase space J _1 (0, T) /Diff vo i(J r o) is symplectically isomorphic to T* SE{2), 
equipped with the left- invariant symplectic form given at the identity by 

Sl r (e):^„(e)-re;Ae; (4.12) 

where w C an is the canonical symplectic form on T* SE{2). 

• The reduced kinetic energy Hamiltonian on T* SE{2) is given by 

H Ied (g,n) = VmT 1 ^, (4.13) 
where M = Mf, + M.f is the full mass matrix. 
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4.2 Reduction by the Euclidian Symmetry Group 

In the previous section, we have eliminated the particle relabeling symmetry by restricting 
to the fluid configurations that had circulation T and no external vorticity. However, after 
reducing by Diff vo i (J-b), the solid- fluid system is still invariant under global translations 
and rotations. In other words, the group SE(2) acts as a symmetry group for the reduced 
equations. This is clear from Theorem 4.6: the group SE(2) acts on T*SE(2) by left 
translations, and both the reduced symplectic form as well as the Hamiltonian are invariant 
under that action. 



The Reduced Poisson Structure. The symplectic structure (4.12) on T*SE(2) is in- 
variant under the left action on SE(2) on itself, and hence induces a Poisson bracket on the 
dual space se(2)*. If the original symplectic structure were the canonical one, the reduced 
Poisson bracket would simply be the (minus) Lie-Poisson bracket (A. 4). However, because 
of the magnetic term in the symplectic structure, additional terms arise in the expression for 
the reduced Poisson bracket. The general form of the magnetic Poisson bracket is described 
in Theorem 3.5. In the case of the rigid body with circulation, the first term in (3.13) is 
the Lie-Poisson bracket on sc(2)*, given by (A. 5). The second term in (3.13) is due to the 
magnetic two-form. The entire Poisson bracket is then given by 

{F,G} B - {F,GU 2) , T ^— — - — — 

To make this bracket more explicit, we may evaluate the Poisson bracket on the coordinate 
functions on se(2)*: 

{P x ,P y } B = -r, {IL,P X } B = -P yi {IL,P y } B = P«. (4.15) 

This emphasizes another advantage of deriving the Poisson brackets through reduction: any 
reduced Poisson bracket automatically satisfies the Jacobi identity, obviating the need for 
any explicit computations. While these computations would have been straightforward in 
this case, this is not always so (see Marsden and Ratiu [1994] for examples). 



(4.14) 



The Equations of Motion. Having established the expression for the reduced Poisson 
bracket on se(2)*, we now turn to the reduced Hamiltonian. Since -ff re d on T*SE(2) in 
(4.13) is written in terms of a left trivialization of T*SE(2), we immediately obtain that 
the reduced Hamiltonian on se(2)* is given by 




(4.16) 



The equations of motion relative to the Poisson bracket (4.14) hence take exactly the form 
given in (2.13) which we repeat here for completeness 

n = (P x v) • b 3 

(4.17) 

p = p x m> 3 + rb 3 x v 

where we have used the fact that 

_ = f2 and _ =V . 

As mentioned in the introduction, these equations occurred first in the work of Chaplygin 
and Lamb, and were given a sound Hamiltonian foundation by Borisov and Mamacv [2006] 
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and Borisov, Kozlov, and Mamaev [2007a]. Moreover, the equations (4.17) are a special case 
of the equations of motion derived by Shashikanth [2005], Borisov, Mamaev, and Ramodanov 
[2007b] and Kanso and Oskouei [2008] for a cylinder of arbitrary shape with circulation 
interacting with point vortices. Our equations are slightly different from the ones derived 
by Chaplygin [1956] and used in Borisov and Mamaev [2006], but can be brought into that 
form by diagonalizing the mass matrix M and rewriting the equations (4.17) in terms of 
velocities rather than the momenta. 

The Kutta-Zhukowski Force. It is worthwhile to point out the geometric significance 
of the equations of motion. For T = 0, the equations (4.17) reduce to the classical Kirchhoff 
equations, which are seen to be Lie-Poisson equations on se(2)* (a more general version 
of this result already appears in Leonard [1997]). When the circulation L is non-zero, an 
additional gyroscopic force appears in the equations of motion, which is traditionally referred 
to as the Kutta-Zhukowski force. In coordinates, this force is given by 

rb 3 xv = r(-^,i4,o), 

and it follows that the force is proportional to P and at right angles to V. From a geometric 
point of view, the Kutta-Zhukowski force arises because of the magnetic terms in the Poisson 
bracket (4.14). Since the magnetic terms are in turn generated by the curvature of the 
Neumann connection, we have hence established the Kutta-Zhukowski force as a curvature- 
related effect. 

A similar effect arises in the dynamics of a charged particle in a magnetic field. Here, 
the particle moves under the influence of the Lorentz force, which is again a gyroscopic 
force whose magnitude is proportional to the charge e of the particle. In the Kaluza-Klein 
approach however, the Lorentz force becomes part of the geometry when the magnetic 
potential is interpreted as a connection whose curvature is precisely the magnetic field 
strength tensor. 

The Sr-potential. We now turn to the computation of the £>f)-potential ip as in defini- 
tion 3.6. For the rigid body with circulation, we will refer to ip : SE(2) — > se(2)* as the 
Bp-potential, in order to emphasize the underlying magnetic form Br- 

Proposition 4.7. The Br-potential tp : SE(2) — > se(2)* is given by 




(4.18) 



Proof. By (3.14), we have that B r (£sE(2), Vse(2)) = • Vse(2) for 6 se(2). The 
infinitesimal generators for the left action of SE(2) on itself are defined by £,sE(2)(g) £, ■ <7> 
and similarly for t]se(2)- Because of the left S'i?(2)-invariance of Br, we have that this is 
equivalent to 

B r (e)(Ad s - 1 e,Ad 9 -x7 7 ) = (d^,r] SE(2) (g)) . (4.19) 
We now write g = (i?e,xo), £ := (SI, V) and i] = (St, V). The vector Ad g -i£ is then given 

by 

Ad (fl , >xo) -i {SI, V) = (n, Rj(V + S2b 3 x x )), 
while the infinitesimal generator t]se(2) is given by 

(O,V) gB(2) CR e ,x ) = ft J, + (V- fijfxo) ■ V. 
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Substituting these expressions into the definition (4.19) then yields 

r r 

- 0(n,v)(^e,x o ) = -O- ||x || 2 + -V T Jx , 
(up to an arbitrary constant, which we set to zero) which is equivalent to (4.18). I 

The Symplectic Leaves in se(2)*. Proposition 3.7 offers a simple prescription for the 
symplectic leaves of the magnetic Poisson structure (4.14). Using the expression (4.18) for 
the £>r-potential, the affine action of SE(2) on se(2)* is given by 

(Rg,K ) • (II, P) = (n - P T Rjlx - j ||x || 2 , R e P + ^Jx ^ , (4.20) 

where (i?g,Xo) € SE(2) and (II, P) G se(2)*. The orbits of this action are paraboloids of 
revolution with symmetry axis the II-axis (see figure 4.1): they are level sets of the Casimir 
function 

$ = n + l||P|| 2 . (4.2i) 

Note how the case with zero circulation is a singular limit of (4.21), and compare this with 
the picture of the standard coadjoint orbits in se(2)* (which can be obtained by setting 
r = in (4.20)), which are cylinders around the II-axis, together with the individual points 
of the II-axis. The trajectories of the system lie in the intersection between the symplectic 
leaves and the level surfaces of the Hamiltonian; see Borisov and Mamaev [2006]. An explicit 
integration by quadratures was obtained by Chaplygin [1956]. 




Figure 4.1: Intersection with the (II, P^-plane of the orbits of the affine action 
(4.20) in se(2)*. 



5 Geodesic Flow on the Oscillator Group 

The dynamics of a rigid body with circulation bears a strong resemblance to the motion of 
a charged particle in a magnetic field. In both cases, the Hamiltonian is the kinetic energy 
and the system is acted upon by a gyroscopic force, either the Kuta-Zhukowksi force, or the 
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Lorentz force. In the case of a magnetic particle, the Lorentz force can be made geometric 
by means of the Kaluza-Klein description. The trajectory of a magnetic particle is then 
a geodesic in a high-dimensional space which is the product of the original configuration 
space M and the group U(l). See Marsden and Ratiu [1994] for an overview and Sternberg 
[1977] for the extension to Yang-Mills fields. 

One can now ask the question whether a similar description exists for the rigid body with 
circulation, such that its motion is geodesic in an appropriate higher-dimensional space. 
Surprisingly, the answer turns out to be positive. The relevant space is the oscillator 
group Osc, which is a central extension of SE(2) by the real line R. Below, we analyze the 
structure of this group and we give an explicit expression for the Lie-Poisson bracket on the 
dual of its Lie algebra. We then show that the equations of motion induced by this bracket 
are closely related to the equations of motion for the rigid body with circulation. 

The use of central extensions in the study of mechanical systems is not new. We mention 
here in particular the work of Ovsicnko and Khcsin [1987] on the KdV equation as a geodesic 
equation on the Virasoro group, and the work of Vizman [2001], in which the structure of 
the geodesic equations on an extension of a Lie group is studied in general. A detailed 
account of the geometry of central extensions, including the oscillator group, can be found 
in Marsden, Misiolek, Ortega, Perlmutter, and Ratiu [2007]. 

The Oscillator Group. The oscillator group Osc is a central extension of SE(2), i.e. 
there is an exact sequence 



-> R -> Osc -> SE(2) -> {e}. 

The multiplication in Osc is determined by the specification of a real-valued SE(2)-two- 
cocycle B : SE{2) x SE{2) -> R: 



Here, Wr2 is the standard symplectic area form on R 2 : 

^R 2 (x,y) := x ■ Jy, 
where J is the symplectic matrix (A. 2). 

The Lie algebra of the oscillator group is denoted as osc and has as its underlying vector 
space the space se(2) x R. The Lie bracket is determined by the Lie algebra two-cocycle 
C : se(2) x se(2) — > R induced by B as follows 



B((R0,x o ), (R^,y )) := w R 2(x ,i? e y ) 



x • IR e x . 



(5.1) 



C(vi,v 2 ) ■ 



dsdt s ,t=o 



(B(g(t),h(s))-B(h(s),g(t))), 



where g(t) and h(s) are smooth curves in SE(2) such that g(0) 
Explicitly, C is given by 



Vi and h(0) = u 2 . 



C((fii s Vi) ) (n 2 ,V 2 )) = 2wR a (Vi s V 2 )=2Vi-JIV 2 . 



(5.2) 



The bracket of the oscillator algebra is then given by 



[(i2i > Vi,a) J (n 2 ,v 2> 6)] = ([(fii J v 1 ) 1 (n 2 ,Va)],c((n 1 ,Vi),(n 3> v 2 )) 
= (o, -r>ijv 2 + o 2 jVi,2Vi ■ jv 2 ). 



Further information about the oscillator group and its algebra can be found in Streater 
[1967] and Marsden, Misiolek, Ortega, Perlmutter, and Ratiu [2007]. 
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The Lie-Poisson Bracket on osc*. We now determine the Lie-Poisson bracket on the 
dual Lie algebra osc* and relate this expression with the Poisson bracket for the rigid body 
with circulation. Note first that as a vector space, osc* is just se(2)* x R, so that an element 
v of osc* can be written as v := (n,p), where it — (II, P) £ se(2)* and p £ R. 



Proposition 5.1. The Lie-Poisson bracket on osc* is given by 

{/,5}o S c*(tt,_p) = {/,g} 5c (2)*(7r) -pC y~ 

\oir Sir 

where C is the sc(2)-two-cocycle given by (5.2). 



(5.3) 



Proof. As usual, it is sufficient to determine the Lie-Poisson bracket on linear functions 
/, g on osc* , for which 

Sf\ c ._ Sf 



for ~ £ osc, 
Sv 



and similar for g. Since v £ osc*, we may write v = (7r,p), where ir £ se(2)* and p £ 
Likewise, we may decompose the variational derivative as 



Sf (Sf Sf 



Sv 



Sir ' dp 



£ se(2) x R. 



The minus Lie-Poisson bracket on osc* is then given by a similar formula as (A. 4), and we 
get 



{/,<?}c 



6£ Sg 
5v ' 8v 



8_l Sg 
Sir ' 5tt 



sf s A 

5tt ' 5tt 



pC 



Sf s A 

5tt ' Sir 

s_i s A 

Stt ' Sit 



so that we obtain the expression (5.3). 



In coordinates, the Poisson structure of proposition 5.1 is given by 



{/,ff}osc*(T,P) = (V^p)/) 5 








~Py 


p. 


°\ 




Py 





-p 







-Px 


P 








V 











0/ 



(V( ff , p )ff). 



(5.4) 



These coordinate expressions also serve to clarify the link between the description of the 
rigid body on the oscillator group and the Euclidian group: each of the subspaces se(2)* is 
a Poisson submanifold of osc* , as in the following proposition. 

Proposition 5.2. Let T £ R. The inclusion Lr : se(2)* <-»• osc* given by tr(i") = (t^T) * s 
a Poisson map. 



Proof. This can easily be seen from the coordinate expressions above: for p = T, the 
Poisson structure (5.3) coincides with (4.14) on the submanifold se(2)*. ■ 
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The Equations of Motion. The reduced kinetic energy Hamiltonian (4.16) on se(2)* 
gives rise to a Hamiltonian H osc * on osc* as follows: 

2 i 2 

H osc *(ir,p) = H(ir) + P - = -n T M" 1 tt + V -. 

The equations of motion for this Hamiltonian and the bracket (5.3) are then given in coor- 
dinates by 

ri = (Px v)-b 3) 

P = P x Qb 3 + pb 3 x V, (5.5) 
p = 0. 

Conclusions. We summarize the conclusions of this section in the following theorem. 
Theorem 5.3. Let V be an element ofM 2 . 

• The Lie-Poisson structure on the oscillator algebra is given by (5.3) and the resulting 
equations of motion by (5.5). 

• The dynamics of a rigid body with circulation T takes place on the Poisson submanifold 
p = r of osc* . 

6 Conclusions and Outlook 

In this paper, we established the non-canonical Hamiltonian structure of Borisov and Ma- 
maev [2006] by geometric means, and we showed how the geometric description gives new 
insight into classical results such as the Kutta-Zhukowski force. We now discuss a number 
of open questions related to this description. 

Cocycles. The problem considered here seems to fit in a general class of systems where 
a remarkable interaction between cocycles and curvature forms is at play; examples include 
the work of Holm and Kupcrshmidt [1988] and Ccndra, Marsdcn, and Ratiu [2004] on the 
analogy between spin glasses and Yang-Mills fluids. Recently, a comprehensive reduction 
theory has been developed for such systems (see Gay-Balmaz and Ratiu [2008] and Gay- 
Balmaz and Ratiu [2009]). It would be interesting to see whether this reduction theory can 
be directly applied to the rigid body with circulation considered here. 

The Oscillator Group and Reduction. It is still not entirely clear why the description 
in terms of the oscillator group, highlighted in proposition (5.2), should exist. Is there 
some way of deriving the oscillator group immediately from the unreduced space Q and 
the particle relabeling symmetry? One way to do this goes back to the work of Kelvin, 
who dealt with circulation by placing a fixed surface in the fluid so that the fluid domain 
becomes simply connected. The flow across the surface then appears as a new variable, and 
the reduced variables are precisely the coordinates on Osc. Kelvin's argument originally 
only held for fluids in a bounded container, but can be extended to unbounded flows as 
well. This approach will be the subject of a future work. 
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Point Vortices and Vortical Structures. In Vankcrschaver, Kanso, and Marsden 
[2009], the dynamics of a rigid body interacting with point vortices was investigated from 
a geometric point of view. It would be of interest to extend the results of this paper to 
the case of a rigid body with circulation moving in a field of point vortices or a general 
distribution of vorticity and, in particular, it is interesting to investigate whether one could 
generalize the oscillator group description to the case of non-zero vorticity. Lastly, we also 
intend to study the geometry of controlled articulated bodies moving in the field of point 
vortices or other vortical structures. This setup has important applications in the theory of 
locomotion; see Kanso, Marsden, Rowley, and Melli-Huber [2005] and Kanso [2009]. 



A The Rigid Body Group 

In this appendix we discuss some elementary aspects of the special Euclidian group SE(2), 
consisting of translations and rotations. This material is standard and can be found, for 
instance, in Marsden and Ratiu [1994] and Arnold, Kozlov, and Neishtadt [1997]. 



The Rigid Body Group SE(2). An element g € SE(2) consists of a rigid rotation and 
a rigid translation: g = (R$,Xo), where Re € 5*0(2) and Xo <E M 2 . The group multiplication 
and inversion in SE{2) are respectively given by 

(i?0,x o ) ■ (i?</,,y ) = (i? e+ ^,x + Reyo) and (R e , x ) _1 = (R-e, -R-exo)- 

It will often be useful to write an element (i?e,Xo) G SE{2) as a matrix 

X 1 °). (A.l) 

Under this identification, the group composition and inversion in SE(2) are given by matrix 
multiplication and inversion. The Lie algebra se(2) of SE(2) is parametrized by Q = (fl, V) T 
and can be identified with R 3 with the following Lie bracket: 

[(fii, V0, (fi a , V 2 )] = (0, -1W 2 + fW 2 ) 
for all f2 12 G K and V 12 € R 2 , where J is the symplectic matrix given by 



1 

-1 



(A.2) 



We will often use the following basis of se(2): 

e n = (1,0, Of, e;c = (0,l,0) T e y = (0, 0, 1) T . (A.3) 

The dual se(2)* of the Lie algebra of SE(2) can also be identified with M 3 by using the 
Euclidian inner product on se(2) = M 3 . We denote a typical element of se(2)* by ir = (II, P) 
with II e 1 and P = (P x ,P y ) £ M. 2 . The notation is chosen to be reminiscent of angular 
and linear momentum in the body frame. The Lie-Poisson bracket 2 on sc(2)* is then given 
by the standard formula: 



{F,GU 2y (,) = -(n,^, S ^ (A.4) 



2 We use the minus Lie-Poisson bracket, which is the relevant choice of Poisson bracket for rigid body 
dynamics. More information about this choice can be found in Marsden and Ratiu [1994]. 
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-Py Px 



for functions F{U, P), G{U, P) on se(2)*, and since 

SF fdF 

we have that the bracket is explicitly given by 

{F,G} se{2r = (V (n ,P)F) T AV (n ,P)G, where A = ( P y 0). (A.5) 

, -P x 

For future reference, we denote by 

e& = (1,0,0), e* x = (0,1,0) = (0,0,1) (A.6) 
the basis of se(2)* = K 3 which is dual to the basis (A. 3) of se(2). 

B The Group of Volume-Preserving Diffeomorphisms 

In this appendix we recall some properties of the group of volume preserving diffeomorphisms 
and its Lie algebra. This material is well-known and can be found for instance in Arnold 
[1966], Ebin and Marsden [1970], and Arnold and Khesin [1998]. 

The fluid domain is a subset of 1R 2 . The reference configuration of the fluid is denoted by 
J-"o and the space taken by the fluid at a generic time t is denoted by T . The fluid domain 
J 7 is a submanifold of ]R 2 and hence it inherits the Euclidian metric from R 2 , which we 
denote by ((v, w)) for all v, w £ TT . The associated Hodge star operator will be denoted by 
* : ^l k {T) — > ft 2 ~ fc (.F), and using * we may write the induced metric on the space of forms 

by 

((a,0)) := / a A*/?. (B.l) 



We also introduce the co-differential S : il k (J-) — > f2* _1 (.F) by the standard prescription 
Sa = * d * a. Lastly, we introduce the index raising and lowering isomorphisms induced by 
the metric, and write them respectively as jj : T*J- — > TJ 7 , where (($(a),v)) = a(w), and 

b = r 1 - 

An incompressible motion of the fluid is a one-parameter (time) family of elements of the 
group Diff vo i(J-o) of volume-preserving diffeomorphisms, which encodes the symmetry of the 
fluid by arbitrary particle relabelings. Formally speaking, DifF vo i(J-o) is a Lie group with as 
its Lie algebra the space X vo i(-?o) of divergence- free vector fields which are tangent to the 
boundary of the body B. The Lie bracket is the negative of the commutator of vector fields. 
The dual space 3£* ol (J"o) can be identified with the space r2 1 (J r o)/dfi (J r o) of one-forms up 
to exact forms with duality pairing 

(X, [a]) = { a(X) dV, (B.2) 



but a better identification can be made. Any equivalence class [a] £ X* ol (To) is uniquely 
determined by the exterior derivative da and by its value on a basis of the first homology of 
J"o- Because of the presence of a single rigid body, the homology is generated by any closed 
curve C around the body, and the value of a on C is simply given by 

a. (B.3) 



c 
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The derivative da corresponds to the vorticity of the system, while T corresponds to the 
circulation. We hence have an isomorphism between X* ol (J-"o) and dl7 1 (J r ) x K given by 

N^(da,r), (B.4) 

where T is given by (B.3). 

We finish by noting that Diffvo^Jo) acts on Xvoi^o) and 3tl ol (Fo) through the adjoint 
and the co-adjoint action, respectively. These actions are both given by push-forward: if 
4> is an element of Diffvo^.Fo), and u and (da,T) are elements of X vo \(J-q) and X* ol (J-"o), 
respectively, then 

Ad (u) = fau, and CoAd (da, T) = (d(0*a), r )- ( B -5) 

The fact that the co-adjoint action leaves the circulation T invariant is related to Kelvin's 
theorem, which says that the circulation around a material loop in the fluid is constant. 
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